Scale-dependent bias from multiple-field inflation 
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We provide a formula for the scaling behaviour of the inflationary bispectrum in the 'squeezed' 
limit where one momentum becomes much smaller than the other two. This determines the scaling 
of the halo bias at low wavenumber and will be an important observable for the next generation of 
galaxy surveys. Our formula allows it to be predicted for the first time for a generic inflationary 
model with multiple light, canonically-normalized scalar fields. 
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For the last twenty years, our information about the 
very early universe has primarily come from the large- 
scale anisotropics of the cosmic microwave background 
(CMB). But over the next decade fresh data will come 
from a new generation of large-scale structure surveys 
such as Euclid and LSST, and complementary small-scale 
CMB observations such as PIXIE [T]. 

These surveys will probe a new suite of observables, re- 
quiring the development of techniques by which they can 
be predicted from models of the early universe. In this 
Letter we focus on an observable which will be a target 
for future surveys — the scale-dependent halo bias, which 
is sensitive to the scaling of the primordial bispectrum in 
the so-called squeezed limit. It is a discriminant of infla- 
tionary models with more than one active field [2J [3] . 

Scale- dependent bias. — We cannot observe the pri- 
mordial density field itself, but only correlations between 
populations of objects which trace its properties. One 
such population are halos of mass M, which form from 
collapse of over-dense regions. The clustering of these 
halos is described by the halo-halo or matter-halo power 
spectra P hh and P , where 8 is the primordial density 
contrast. At lowest order the relationship between 8 and 
the halo number density is linear P Sh (k) ~ b\ J Pf^{k) and 
therefore P hll (fc) blP R s {k). The superscript 'i?' de- 
notes smoothing over a lengthscale R sufficiently large to 
enclose the halo mass M. 

Corrections to &l come from higher-order correlations 
of 8. They can be expressed using the n-point functions 
of the curvature perturbation which is typically used 
to characterize predictions of inflationary models. The 
leading effect comes from the three-point function and 
induces a fc-dependent shift &l — > bh + Sb(k), where [H[S] 



Sb(k) 
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Here, z is the redshift and Q 2 



k 2 



critical density for spherical collapse is 8 C 
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a\j is the variance of the tracing population of mass-M 
halos [6]. The kernel M R (q, z) is related to the linearized 
Poisson equation connecting the density contrast (5(k, z) 
and C(k). It can be written as 



M(k, z) 



2 k 2 T(k) g(z) 
5flffi^l + z' 



(2) 



The transfer function T(k) is normalized to unity as k — ¥ 
0, and g(z) accounts for suppression of growth during 
A-domination. The present-day matter density is f2 mj o- 
To smooth <5(k, z) on the scale R we introduce a window 
function W R {k) and define M R (k,z) = W R (k)M[k,z). 
Finally, the spectrum and bispectrum of the curvature 
perturbation satisfy 



(C(k 1 )C(k 2 )) = (27T) 3 (5(k 1 
(C(k 1 )C(k 2 )C(k 3 )) = (27T) 3 <5(k 1 



k 2 )P c (3) 
k 2 +k 3 )B c . (4) 



Scaling of < 



-For k smaller than the wavenumber k t 



eq 

corresponding to the horizon scale at matter-radiation 
equality, we have T(k) ~ 1 and M(k, z) proportional to 
k 2 . For k < k cq this sets the dominant scaling of 8b. It is 
the complete result whenever the integral over Bq does 
not depend on k, as for the case of the local model with 
constant amplitude /nl [H E] • 

In this Letter we focus on the possibility that the in- 
tegral in ([T]) also scales nontrivially with k. The precise 
scaling depends on details of the underlying inflation- 
ary model, and may constitute a useful observable. Be- 
cause M(q) grows steeply with q for q < fc cq , and the 
window function restricts the integral to wavenumbers 
q < R , the dominant contribution will come from the 



region k, 



eq 



< 
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If k -C k cq then this is the 



'squeezed' region, where one of the momenta in the bis- 
pectrum is much less than the other two. Therefore, for 
k <C fc eq , the halo bias is a probe of the squeezed limit of 
the bispectrum. 

To study this limit we denote the 'squeezed' momen- 
tum k L . The other two momenta are taken to be of order 
k s , with fc L <C k s . We will see below that it has not 
yet been possible to obtain reliable predictions for the 
bispectrum in this limit, because the presence of multi- 
ple large hierarchies causes perturbation theory to break 
down. But if k h is not too much smaller than k s , the 
bispectrum can be approximated by a power law with 
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spectral index n sq , 

£<(fcs) 
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where B ( -(k s ) is an overall amplitude depending weakly 
on k s and k t = fcj is the perimeter of the momentum 
triangle ki+k2+k3 = 0. The spectral index n sq measures 
the response of to a change of k h and may also depend 
weakly on k s . In what follows, our aim is to provide a 
precise formula for the spectral index n sq and relate it to 
the scaling of Sb. 

The integral over q in (fil has been obtained for a va- 
riety of simple templates [3l [7] • Here we focus on the 
scale dependence inherited from the squeezed limit of the 
bispectrum. The amplitude and scale-dependence of the 
spectrum and bispectrum can be expressed in terms of 
an arbitrary reference scale k* . We write 



P c (k) 



2fc 3 



2fc 3 



n,(k,)-l 
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where n s — 1 = dlnT^/dln/c is the spectral index of P^, 
and B ( -(k s ) = £^(/c*)(£; s /fc*) Q . We assume n s and n sq to 
be approximately constant over the range fc eq < q < R . 
Because this is comparatively small, with Rk eq ~ 0.1 
on cluster scales, we expect this to be reasonable. We 
conclude that Sb has approximate scaling behaviour 



5b « A Nh x M R (k,z)~ 



-(n.-l) 



(7) 



where the amplitude Anl depends on details of the 
model, the scale R and the background cosmology. In 
general, it may be difficult to calculate. For the purposes 
of this paper we do not require a precise estimate, be- 
cause the observable is the scaling behaviour of Sb rather 
than its amplitude. This is a strength of the approach. 

We define the combination nst, = n sq — n s + 1 as the 
spectral index of the bias. The main result of this Letter 
is a prescription to compute it for any inflationary model. 
Constraints on rigb have been obtained from present-day 
data [2], although the results show sensitivity to priors. 
In particular, the first Planck results impose strong con- 
straints on the value of /nl local [5] and therefore on the 
range of models that can show an observable ngb- Im- 
provements are expected from a future Euclid- or LSST- 
like survey [9]. 

The fc t -dependence of the bispectrum Bq was first stud- 
ied by Chen [10] . Formulae for a multiple-field model 
were given by Byrnes et al. [TTJ [T2]. Our analysis dif- 
fers because of its focus on scaling in the squeezed limit, 
rather than scaling of the amplitude for nearly equilat- 
eral momenta. Later, Tzavara & van Tent studied the 
/^-dependence and the scaling behavior of ^ in a two- 
field model using a Green's function formalism 1 1 3 j . In 
this Letter, we give expressions using the separate uni- 
verse formalism, in both its 'variational' and 'transport' 



versions, which are valid for any number of scalar fields. 
These expressions are only valid when super-horizon evo- 
lution of curvature perturbations has ceased, in other 
words, when the adiabatic limit has been reached. Oth- 
erwise, if isocurvature perturbations persist at the end of 
inflation, it is necessary to augment the model by pro- 
viding a prescription for their subsequent evolution. Our 
formulae apply when it is sufficient to evaluate the curva- 
ture perturbation at or before the end of the inflationary 
phase. 

The scalar spectral index. — Sasaki & Stewart used 
the separate universe method to give an expression for 
n s [14] . More recently an alternative (but equivalent) 
prescription was given in Ref. [15] . which we follow here 
[16] . According to the separate universe picture, the 
power spectrum evaluated at time t can be written as 



V c {k)\ t = N a NpE aP (k)\ t0 , 



(8) 



where N a = dN(t, to)/<90 Q (to), and N measures the 
number of efolds elapsed from a spatially flat slice at 
time to to a uniform density hypersurface at later time 
t; indices a, (3, . . . , label the different species of scalar 
field; and E Qj g defines the two-point function of scalar 
field fluctuations, 



(^ a (k!)^(k 2 )) t0 = (2 7 r) 3 «5(k 1 + k 2 ) 



2fc 3 



(9) 



with fc = |ki| = | ki2 1 - The time to can be chosen to 
coincide with the horizon crossing time of k, which we 
will denote with a subscript l k\ 

The spectral index can be obtained from ^ provided 
we know to next-order in slow-roll, which was first 
obtained by Nakamura & Stewart [T7]. At this order and 
for light, canonically-normalized fields, we find 



Hi I S, 



a/3 



2j 



a/3 



2< /3 l n (-fc*r) 



(. 10 ) 

where '*' denotes evaluation at the horizon-crossing time 
for an arbitrary scale the conformal time r satisfies 
dt = a(t)dr; r a /3 is a constant; u a $ = —m a p/3H 2 is a 
rescaled mass matrix for the fluctuations, equivalent to 
the expansion tensor of the inflationary flow field [18] ; 
M a p = eS a/ 3 + Uap, and e = —H/H 2 . 

The fc-dependence can be extracted from the coefficient 
of lnfe/fc*, provided we know how to choose fc*. To do 
so, note that the ln(— fc*r) term is the lowest power in 
a series expansion which describes the time dependence 
of the fluctuations. Since we wish to estimate the scale 
dependence at a fixed time to we can choose the arbitrary 
scale so that fe*r = —1, making all powers of ln(— fc^r) 
negligible. With this choice the evaluation point of H, 
Tap, Uap and M a p becomes coincident with to and the 
spectral index simplifies to 
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t N a NpM aP \ k 
' N X N X 



(11) 
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Squeezed bispectrum. — A similar approach can be eral single-field models [THj, and were recently obtained 

used to obtain the spectral index n sq in ([5| . We require for multiple canonically- normalized fields in Ref. [20) . In 

a next-order expression for the bispectrum, in analogy the squeezed limit they can be written 
with (10). The requisite corrections are known for gen- 



( (k L )<%(k s )^ 7 (k' s )) = (27r) 3 <5(k L + k s + 14)5, 



a/37 



fc L ->0 



(27r) 3 £(k s + 14 
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(12) 



Corrections of order ~ k^/kg may be relevant in the 'not-so squeezed' limit where the hierarchy k L /k s is modest |21j . 
but in the present case we expect their effect to be negligible. The coefficient b a \p^ is symmetric under interchange of 
the indices /3 and 7, but not necessarily under other permutations. In what follows it is helpful to break it into terms 
of lowest-order and next-order in slow-roll, which we label 'lo' and 'nlo'. Using the results of Ref. [20], we find 



LO _ 


tt4 H>u x 


NLO 





LLO 

5 A|/3 7 



«■ A ^ A ) ln(-fc»r) - 2M* x b^ In 



1LO 

?A|/3 7 - 



2fcL 



(13a) 
(13b) 



The omitted terms are not logarithmically enhanced, 
or are proportional to ln(A; s /fe*) or ln(fct/fc*). In (13a) 
and (13b) we have chosen fc* ~ k t , making these contri- 
butions negligible Therefore, if we intend to estimate this 
expression at some fixed time, only ln(fc L //c*) can gener- 
ate dangerously large contributions. In analogy with the 
ln(— fc*r) term in (101, this is the lowest term in a series 
expansion, and, if uncontrolled, the series will diverge 
when fc L /fc* R3 k h /k s — > 0. It is for this reason that there 
is no analytic formula for the squeezed limit of the bis- 



pectrum in a multiple-field model. 

Ref. [20] developed a method based on the dynamical 
renormalization group to control the series of ln(— fc*r) 
terms in the limit t — > 0. It is not yet clear whether 
a similar scheme could be used to deal with the series 
expansion of hi(k h /k s ). For this reason our results are 
untrustworthy for fc L <C k s . Our approach may be reliable 
up to a hierarchy | ln(fc L /fc s )| of order a few. 

Using the separate universe formula, the bispectrum 
-Bj can be written (for an arbitrary momentum triangle) 



B ( (ki,k2,k 3 )\t = N a NpN 7 B a p 7 (ki,k 2 ,k 3 )\ to 

"S aj (ki) E/3s(k 2 ) S Q7 (fci) T>ps(k 3 ) 



+ N a pN^N S 



E Q7 (fc 2 ) £/3s(fc 3 ) 



"2 



1,3 

"•3 



1,3 

lx 2 



h 3 



(14) 



where N a p = d 2 N{t, to)/d(f) a (to)d(f)p(to). The squeezed spectral index, n sq , can be computed using this expression. 
We take to to be the horizon crossing time for k t and denote evaluation at this time with a subscript 'kt'. Choosing 
fe* ~ kt, the term proportional to ln(— fc*r) in (13b I vanishes and the scale dependence of B a p^ can be read off. In 
combination with (11), we obtain 



( NaNpNry (M a A& A ° / 3 7 ) I kt + N af) N y Ns(M a xZx'y + M 7A S aA ) | fct S^| fet 
! N x N^N v b^J kt + 2N Xll N u N^ Xu V^)\ kt 



(15) 



r 



Whenever the bispectrum in a multiple-field inflation- 
ary model is large enough to be observable, the 'nonlin- 
ear' terms in the second line of (fl4| dominate the 'in- 
trinsic' term in the first line [351 123j . In this limit, since 
^ap\k t ~ B\ $ a/3, we can write n sq in the simpler form 



(16) 



and under the same assumptions we conclude that 

U5b ~~ 2 N X ,N X N, +2 N X N X ■ (17) 
This constitutes one of the principal results of this Letter. 



Transport equations. — Eqs. (11) and (15) are framed 
in terms of the 'variational' formulation of the separate 
universe method, which involves 'variational' derivatives 
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such as N a . A successful numerical implementation re- 
quires sufficient resolution to extract the variation of N 
against a background of numerical noise accumulated 
during the integration from to to t. An alternative ap- 
proach to solving (l8|) and ( 14 ) consists in integrating 



the correlators of field fluctuations using Jacobi evolu- 
tion equations [TS]. Mathematically, both methods are 
equivalent, but the Jacobi approach is simpler because it 
is more tolerant to numerical noise. For details on this 
formalism we refer to Refs. [181 . 124H28] . 

To write a Jacobi-type equation which defines the time 
evolution of the coefficient matrix 6 Q | ( g 7 in ( 12 ), we define 
u a py = du a /3 1 d(f>-y . Using the prescription of Ref. |20j , it 
follows that 



db, 



dlnr 



(18) 

Note that we should regard one of the two-point functions 
E a/ 3 to be evaluated at k h and the other at fc s . A suitable 
initial condition at horizon-crossing is provided by the 
lowest-order value (13a), after setting ~ k t . 



To compute the scale dependence in the squeezed limit, 
we adjust fc L while keeping k s fixed. Defining n a \p^ = 
db a \pj/dhik h , it follows that 

dn a\P~f _ 

dhir ~ aX X \P~I U P* n a\\y + u iX n a\f3X /^g% 



where n a p = dS Q( g/dlnfc. The initial condition is given 
by n Q | / 3 7 = — 2M a xb T £tp r y near horizon-exit for k t . 
Likewise and n a p satisfy [TS] 



d£ 



a0 



(20a) 
(20b) 

with respective initial conditions Y, a p = H^5 a p and 
n Q( g = —M a \E\p — Mpx'EaX near the horizon-crossing 



dN 
dn a p 
dN 



u a xnxp + Ufsxn a x, 



time for k t , supplied by (10) 



To relate n a p and n a \pj to the spectral indices n s and 
n sq requires the gauge transformation to £. We label 
these coefficients iV* because they are obtained by com- 
puting a perturbation in the e-folding number N, but the 
label t emphasizes that they involve only local quantities 
at t. Explicit expressions for these gauge transformations 
are tabulated in Anderson et al. [25] . The spectral index 
of the bias can be written 



nsb 



N^Np^n^ + 2NqNpN^ s n aj Y,ps NjN^ 

(21) 

Focusing on multiple-field models with a bispectrum 
large enough to be detectable, in analogy with Eq. (16), 
this expression reduces to 



(22) 




f NL =0.007 



25 30 35 40 45 50 



(a) Double quadratic model 




(b) 4 /axion model 



FIG. 1: Evolution of ngb some cfolds before the end of inflation in 
two different models. The values of /nl shown are esti- 
mated at the end of inflation. 



As an example, Fig. [T] depicts the evolution of nst in 
double quadratic inflation [29-35] and in the 4 -plus- 
axion model studied by Elliston et al. [33J . To generate a 
significant nsb, the bispectrum must scale differently to 
the power spectrum. If |/nl| is small this is typical, as 
in Fig. la Alternatively, if |/nl| 



is 0(1) or larger the 
scaling is more similar, as in Fig. |lb[ Where the spectral 
index of the bias is large we find that it is often negative, 
giving the squeezed limit a red tilt. This is the opposite of 
quasi-single-field inflation (QSFI), where ngb > 0. This 
can be understood heuristically. QSFI contains mas- 
sive modes, which mediate finite-range forces and tend 
to soften long-range correlations. In comparison, models 
with multiple active fields tend to add tachyons to the 
spectrum which mediate long-range forces and therefore 
enhance correlations. 

Conclusions. — The halo bias is known to inherit a scale 
dependence from the underlying inflationary model. It is 
therefore an important observable which is likely to be 
measured by future surveys of the galaxy distribution. 
In this Letter we focus on its scaling, for which we do 
not require detailed knowledge of the amplitude Anl- 
We provide a formula to compute n$b for any model of 
multiple-field inflation. To obtain a significant effect, the 
squeezed limit (characterized by the spectral index n sq ) 
should scale differently to the power spectrum (charac- 
terized by n s — 1). This can be measured most easily 
when nsb < 0, making k~ 2+nsb diverge more strongly in 
the limit k — > 0. Therefore, a slightly red-tilted bias will 
be easier to constrain than a blue-tilted one. 

The principal drawback of our method is the use 
of a spectral index to parametrize the behavior of the 
squeezed bispectrum. For a sufficiently large variation 
of fc L the bispectrum may have a shape which cannot 
be approximated by a power law. Nevertheless, our for- 
mula should give a good approximation provided the hi- 
erarchy ln(fc L //c s ) does not exceed ~ a few. Near-future 
surveys such as DES may probe a hierarchy of order 
ln(fc L /fc s ) ~ —2, but surveys arriving in the medium- to 
long-term, such as Euclid, may probe ln(/c L /fc s ) ~ —8. 
Small-scale observables such as ^-distortions may even 
probe ln(A; L /fc s ) ~ —19. For such cases a more precise 
description of the bispectrum will be required, perhaps 
using numerical methods. 
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